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The elastic, thermodynamic, and electronic properties of fluorite RuC>2 under high pressure are investigated by 
plane-wave pseudopotential density functional theory. The optimized lattice parameters, elastic constants, bulk 
modulus, and shear modulus are consistent with other theoretical values. The phase transition from modified 
fluorite-typeto fluorite is 88 GPa (by localized density approximation, LDA) or 1 15.5 GPa (by generalized gradient 
approximation, GGA). The Young's modulus and Lame's coefficients are also studied under high pressure. The 
structure turned out to be stable for the pressure up to 120 GPa by calculating elastic constants. In addition, 
the thermodynamic properties, including the Debye temperature, heat capacity, thermal expansion coefficient, 
Gruneisen parameter, and Poisson's ratio, are investigated. A small band gap is found in the electronic structure 
of fluorite RuC>2 and the bandwidth increases with the pressure. Also, the present mechanical and electronic 
properties demonstrate that the bonding nature is a combination of covalent, ionic, and metallic contributions. 
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1. Introduction 

Attractions to study R11O2 are due to its fundamental properties and potential superhard character- 
istics . An extensive search for new superhard materials has been undertaken during recent years 
and a new class of hard materials has been suggested: the transition-metal dioxides containing heavy 
elements. Typically, the bulk modulus of modified fluorite (pyrite phase, Pa3) RUO2 was found to be 
399 GPa 1 3], which is the highest value except for diamond (442 GPa) |4]. Ru0 2 has the rutile (P4 2 /mnm) 
structure under usual conditions |3], and can be transformed to an orthorhombic (CaCl2-type, Pnnm) 
structure at about 6 GPa |5] or 11.8 GPa |6] and to a pyrite structure at about 12 GPa |2, 5]. Moreover, the 
theory indicates the Pa3 structure can be transformed to a fluorite (Fm3m) structure at about 89 GPa or 
97 GPa d. 

Recently, elastic properties focusing on Pa3 phase of RUO2 have been investigated systematically jULll. 
ri I — 1 r — I r - l j — l 
Electronic structures 1 8-11] and optical properties H2hH4Jl of rutile and orthorhombic 1 14] RUO2 have 

been extensively studied. A full-potential linear muffin-tin orbital calculation on the electronic structure 

and bulk modulus of fluorite RUO2 has also been performed flal . The hardness and elasticity in cubic 

RuG"2 and Raman scattering of the rutile-to-CaCl2 phase transition have been probed experimentally 116]. 
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Previous investigations on fluorite Ru02 are not complete and some problems remain unresolved. 
Many properties, such as the hardness, stabilization, elastic and thermodynamic properties etc under 
high pressure, are still unknown. To reveal the superhard characteristics appropriately, a detailed theo- 
retical description of the elastic and electronic properties is necessary. 

2. Theoretical approaches 

In this work, all the calculations have been performed with CASTEP In the electronic structure 

calculations, we have used the non-local ultrasoft pseudopotential 1 19], together with the revised Perdew- 
Burke-Ernzerhof (RPBE) generalized gradient approximation (GGA) exchange-correlation function 1 20] . 
Considering the computational cost, a plane-wave basis set with an energy cut-off of 600.0 eV 0] has 
been applied, and the 12 x 12 x 12 Monkhorst-Pack mesh has been used for the Brillouin-zone (BZ) /c-point 
sampling. Pseudo atomic calculations have been performed for Ru (4s 2 4p 6 4d 7 5s 1 ) and (2s 2 2p 4 ), where 
the self-consistent convergence of the total energy is at 5.0 x 10~ 7 eV/atom. 

3. Results and discussion 

In the equilibrium geometry calculations of fluorite RUO2, both the GGA and the LDA methods have 
been used. The bulk modulus (Bo) and its first-order pressure derivative (B' Q ) by Murnaghan 1 21] equation 
of state are listed in tablejl] 

Table 1. The calculated lattice constants a (A), phase transition pressure Pt (GPa) and elastic constants 
en, C44, c\2 (GPa) by LDA and GGA methods at GPa and K. The bulk modulus B and shear modulus G 
are calculated by the elastic constants. The bulk modulus Bo (GPa) and its first-order pressure derivatives 
B' Q are fitted by the Murnaghan equation of state. 
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The LDA/CAPZ (embedded in CASTEP) calculations are performed using the same parameter input with GGA. 



Figure[T]suggests that a significant stiffer compressibility of 82.92% is obtained using LDA at 100 GPa 
as compared to GGA (80.27%). The volume compressibility (about 80%) is nearly the same as the ultra- 
stiff cubic Ti0 2 at the same pressure |24]. The bulk modulus of the fluorite Ru0 2 is 3S3 GPa (LDA) or 
287 GPa (GGA), which is slightly smaller than that of the Ti0 2 (282 GPa 395 GPa). As the pressure is in- 
creased to 100 GPa, the approximation value of the normalized volume of diamond, c-BN, OSO2, and OsC, 
is 85%, 83%, 80%, and 84%, respectively 1 25-27], which is slightly larger than the current calculations. As 
a comparison, figure [1] shows the Ru-0 bond length contraction with increasing the pressure. Similarly, 
a smaller Ru-0 bond length contraction is obtained using LDA. In a word, the current Ru-0 bond length 
contraction, using either LDA or GGA, is much smaller than that in volume contraction. 

The calculated elastic constants of fluorite Ru0 2 under different pressures are listed in table [2] Ac- 
cording to the generalized elastic stability criteria 1 28] (en - ci 2 > 0, c\\ + 2ci 2 > 0, C44 > 0) for cubic 
crystals, we can demonstrate that fluorite structure is elastic stable under 120 GPa. The larger Cu and 
smaller c\2 indicate the inter-atomic bonding along the c-axis stronger than that along the a-axis, consis- 
tent with the case of the larger bulk modulus B and the smaller shear modulus G over a wide pressure 
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Figure 1. Variation of the lattice volume (denoted by square and circle symbols) and the Ru-0 bond length 
(denoted by triangle symbols) with pressure. They are normalized by XI Xq, where X and Xq are the 
lattice volume or Ru-0 bond length at any given pressure and zero pressure at zero temperature. 

range. Compared to c\2 and C44, C\\ varies largely by changing the pressure, meaning that it is more 
difficult to obtain the same strain from the longitudinal direction than from the transverse direction. 

An estimate of the zero-temperature transition pressure between the Pa3 and Fm3m structures may 
be obtained from the usual condition of equal enthalpies, i.e., the pressure P, at which enthalpy H - 
E + PV of both phases is the same. Our calculated Pa3^Fm3m phase transition pressure is 115.5 GPa by 
GGA and 88 GPa by LDA, as shown in figure[2] in accordance with the theoretical values of 89 GPa |2] and 
97 GPa d but is greater than the predicted 65 GPa flllijl. 

Using the calculated elastic constants at K and GPa, we obtain the bulk moduli B of fluorite RuC>2, 
with the values of 287 GPa (GGA) and 353 GPa (LDA), respectively, which is smaller than that of diamond 
(442 GPa 1 27]), although both of them have comparable compressibility at 100 GPa. The shear constant C44 
is 170.20 GPa (GGA), which is consistent with previous theoretical calculations 152 GPa 12], 140 GPa fill . 
147 GPa 11611 . and experimental measurement 144 GPa 1 16], but is larger than the other theoretical value 
62 GPa l23ll . In general, the shear modulus of cubic materials is slightly lower than the value of C44 fltil . 
whereas our calculation indicates the opposite case. 

There have been proposals that the shear modulus may be a better index of hardness | 29]. Our calcu- 



lable 2. The calculated (by GGA method) elastic constants c\\, c 44 , cyi (GPa), heat capacity Cy 
(J- mol _1 -K _1 ), Debye temperature O (K), Gruneisen parameter 7, thermal expansion coefficient a 
(10 -5 K -1 ) and Poisson's ratio a over a wide pressure range at zero temperature. 
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lated values of G, 225.98 GPa by LDA and 198.16 GPa by GGA, are much greater than the theoretical values 
of 133 GPa |2] and 62 GPa 1 23] . Even so, the current results are still incomparable with those of diamond 
and c-BN, e.g., the recently calculated values for diamond are 550 GPa 545 GPa (H, 518 GPa filll . 
and 403 GPa for c-BN 1 30]. However, the current results are comparable to those of Os02 1 31] with val- 
ues of 250 GPa (LDA) and 223 GPa (GGA) CI 15] suggests Os0 2 as a better candidate for a hard material 
since their calculations confirmed that the bulk modulus of OSO2, 452 GPa, is only smaller than that of 
diamond). Overall, it is reasonable to suggest that the fluorite R.UO2 is a potential ultra-incompressible 
material, consistent with the suggest superhard material from (2J. 




Pressure(GPa) Pressure(GPa) 



Figure 2. Enthalpy as a function of pressures for 
the Pa3 and Fm3m phases of RuC>2, (a) is the LDA 
results and (b) is GGA ones. 



Figure 3. Pressure dependences of mechanical 
quantities (by GGA method) under different pres- 
sures, X represents Bulk modulus B, Shear mod- 
ulus G, Young's modulus Y, Lame's coefficients k. 



In order to better understand the pressure responses of the mechanical behavior, we have studied the 
bulk modulus B, shear modulus G, Young's modulus Y, and Lame's coefficients A by increasing the pres- 
sure to 150 GPa. The Young's modulus Y and Lame's coefficients A are also essential for understanding the 
macroscopic mechanical properties of solids and for designing hard materials. Figure [3] shows the most 
significant pressure dependence of B and the least significant pressure dependence of G. In contrast to B, 
the Lame's coefficients A increase slowly with pressure. Compared to Lame's coefficients A, the Young's 
modulus Y behaves much slower with the increase of pressure. At zero pressure, the relative magnitude 
of the four mechanical parameters in descending order is: Y > B > G > A. However, the Lame's coeffi- 
cient A is larger than G above nearly 30 GPa, and the B is larger than Y above nearly 150 GPa. The high Y 
and B, particularly at high pressures, also suggest that fluorite Ru0 2 is a potential ultrahard material. 

The value of the Poisson's ratio for covalent materials is small (a = 0.1), whereas for ionic materials, 
a typical value of a is 0.25 [32]. In our cases, the value of a for RuC>2 varies from about 0.2124 to 0.3227, 
as shown in table [2] indicating a higher ionic and weaker covalent contribution to intra-atomic bonding. 
Besides, the typical relation between bulk and shear modulus is, respectively, G ~ LIB and G « 0.6B 
for covalent and ionic materials. In our cases, the calculated values of GIB are in the range of 0.7396 at 
GPa to 0.3656 at 150 GPa, indicating that the ionic bonding is dominant for fluorite R.UO2. To evaluate 
the material ductility or brittleness, Pugh et al. introduced the BIG ratio liill : the material is brittle if 
the ratio is less than the critical value 1.75. Therefore, fluorite RUO2 is brittle under ambient conditions 
since the Bl G is only 1.35. However, the brittleness decreases (or ductility increases) when the pressure 
is increased and the BIG ratio rises to 2.4 when the pressure is up to 120 GPa. 

The dependences of the Debye temperature 0, heat capacity CV, Griineisen parameter 7, thermal ex- 
pansion coefficient a, and Poisson's ratio a on the pressure are calculated. As shown in table[2] when the 
temperature keeps constant (T — K), and a increase with increasing the pressure, whereas CV, 7, and 
a decrease. The five thermodynamic parameters show different pressure dependences within the range 
of -r 120 GPa. It is obvious that the thermal expansion coefficient a declines most significantly, corre- 
sponding to an 85% compression. The heat capacity CV and Griineisen parameter 7, however, correspond 
to smaller compressions with 60% and 25%, respectively. The other two parameters, Debye temperature 
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and Poisson's ratio a, increase with the pressure with 32% and 52%, respectively. Moreover, all the five 
parameters have shown decreased dependences with increasing the pressure, indicating anharmonicity 
of the vibration. Therefore, it is necessary to further investigate the electronic energy band structure and 
density of states (DOS) to better understand the physical properties. Accordingly, we have made a sys- 
tematic investigation of the fluorite RUO2 at different pressures (0, 30, 60, 90 GPa) under K, as shown in 
figures H [5] and [6] 




w 



W K0 40 80 
Pressure (GPa) 



Figure 4. Energy band structure (by GGA method) along the high symmetry points in the Brillouin zone 
at the pressures of 0, 90 GPa is shown in (a) and the band gap Eg as a function with the applied pressures 
is shown in (b). 




Energy (eV) 

Figure 5. (Color online) Partial density of states (by GGA method) of O and Ru states under P = 0, 30, 60, 
90 GPa. 

Figure|4]presents the pressure-induced energy level shift towards higher and lower regions. We can 
see that the applied pressure has a larger effect on the energy levels far away from the Fermi level than 
those in the vicinity of the Fermi level, indicating a stronger effect on the core level than on the valence 
level. From the energy band structure, we find that the top of the valence band occurs at W point and 
the bottom of the conduction band occurs at L point (slightly lower than X point by a value of 0.02 eV), 
implying that there exists an indirect gap with width of 0.5175 eV in fluorite RUO2. The calculated band 
gap at zero pressure is consistent with the other theoretical value (0.5 eV) |2], but is much smaller than 
those of diamond (4.15 eV) and c-BN (4.49 eV) [34]. Moreover, it is found that the band gap increases 
almost linearly with the pressure. 
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Figure 6. (Color online) Total density of states (by GGA method) under P ■ 
Fermi level. 



0, 30, 60, 90 GPa near the 



In figures [5] and [6] we plot the calculated DOS by GGA, where the Fermi energy is taken to be zero. 
From figure[5](a), the O 2s band center is at -18.5 eV, which is consistent with those in rutile 1 10, 14] and 
CaCl2-type 1 14] Ru02 . The valence band width is about 4 eV, which is much larger than those in rutile (lj|l 
(2.5 eV) and CaCl 2 -type Q Ru0 2 (1.4 eV). In figure [5] (b), the calculated valence band width of 2p is 
about 8 eV, which is slightly larger than that in rutile Ru0 2 of 5.9 eV (using an extended linear augmented 
plane wave potential) and 6.8 eV (using linear-muffin-tin-orbital potential) 1 10], but the present calcula- 
tion is consistent with that in rutile 1 14] and CaCl 2 -type 1 14] Ru0 2 with the same values of about 8.1 eV. 
Furthermore, the calculated conduction band width is 4 eV, which is far smaller than that of valence 
conduction. 

The Ru s semi-core band, centered at -73.3 eV, displays larger relative intensity (3.6) and smaller 
width (1.2 eV) with respect to those in the conduction band with smaller relative intensity (1.9) and larger 
width (2.2 eV). The other Ru s electrons distribute mainly in the ranges of -21 eVn- -18 eV, -8.6 eV + 
-1.6 eV, -1.5 eV + 0.5 eV, have formed very weak peak with intensities less than 0.1, and thus could be 
ignored as compared with the peak far away from the Fermi level. In figure[5](d), the Ru p state locates at 
-43.5 eV in the valence band with width of 1.3 eV and distributes in the energy range of 6.3h-27.5 eV in the 
conduction band (corresponding to two sharp peaks centered at around 10 eV and 16 eV, respectively). 
The relative intensity of Ru p state in the valence band is far greater than that in conduction band, and 
the DOS distributed in the energy range -20+5 eV is ignored due to their subtle relative intensities (below 
0.12). The Ru 4d state, shown in figure[5](e), is distributed at -18.8 eV and around the Fermi level with 
widths of 2.3 eV and 12.0 eV, respectively. Moreover, the relative intensity of the inner valence band is 
considerably smaller than that in the outer valence band and in the conduction band. 

The four sharp peaks of the total DOS within -7.5 5 eV originate from the strong hybridization be- 
tween Ru 4d and O 2p, as seen in figure[6] The complete overlap of the Ru 4d and O 2p states from -8 eV 
to 4 eV indicates a strong covalent interaction between them, whereas the nonzero DOS value at Fermi 
level indicates a moderate metallic feature in its bonding state. Although there is a large hybridization 
between Ru 4d and O 2p states, the charge transfer from Ru to O is possible in the present case. By ana- 
lyzing the Mulliken population results, it is found that the charge transfer from Ru to O is as numerous 
as about 1.01 electrons. Therefore, the bonding behavior between Ru-0 has ionic contributions owing to 
charge transfer. In a word, the bonding behavior between the Ru-0 is a combination of covalent, metallic 
and ionic contributions. 

To emphasize the pressure dependence of the DOS, we have investigated the DOS under different 
pressures (30, 60, 90 GPa), as shown in figures [5] and [6] It is clearly seen that the applied pressure causes 
the energy levels shifting towards both sides of the Fermi level and thus the energy band is broadened. 
Under a higher pressure, the energy level shift is decreased, implying the strong repulsion among the 
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core electrons. Meanwhile, in general, the relative intensity decreases with increasing the pressure. Ac- 
cordingly, the relative shift in the lower energy space is always larger than those in the higher energy 
space, implying different bonding strength. The change of the DOS can be attributed to the charge trans- 
fer during lattice distortion. With the increase of pressure, a higher overlap of the wavefunction results 
in a stronger derealization of electrons. Electrons transfer from the majority to minority spin band and 
form broader bands. The center changes and electrons become more localized when lattice distortion 
changes from negative to positive. The majority and minority bands move with respect to the Fermi 
level, which affects both the spin polarization ratio and magnetic properties. The current investigations 
reveal that the relative intensity of s and p states decrease with the pressure both in the valence and 
conduction bands. However, the relative intensity of Ru s state keeps almost unchanged in the semi-core 
band and the main peak in the conduction band has been split into two peaks with the pressure. By 
analyzing Ru p state, we find that the relative intensity decreases slightly with the pressure in the whole 
valence band and in the higher-energy range from 12.5 eV to 25 eV, but increases with the pressure within 
7.5 eVn-12.5 eV. Interestingly, there is observed an increase of the relative intensity of Ru d state with the 
pressure in the deeper-lying valence band, whereas the relative intensity of the other Ru d state decreases 
with the pressure, presenting opposite variation tendencies. The different intensity variation trends have 
unambiguously demonstrated that the applied pressure has induced various charge transfer tendencies. 

4. Conclusions 

The current investigations revealed that the fluorite RUO2 is a potential ultrahard material. The elastic 
stability criteria show that the fluorite RUO2 is elastic stable up to 120 GPa. The calculated Poisson's ratio 
and Debye temperature increase monotonously with the pressure. However, the heat capacity, Griineisen 
parameter, and thermal expansion coefficient decrease with the pressure. An analysis based on Poisson's 
ratio, GIB, and DOS reveals that the bonding nature in RUO2 is a combination of ionic, covalent and 
metallic contributions, which contributes to the hardness and fundamental properties. The energy band 
investigations found an energy gap between the top of the valence band and the bottom of the conduction 
band, and the gap seems to increase monotonously with the pressure. Moreover, the different intensity 
variation trends of DOS have unambiguously demonstrated that the applied pressure has caused various 
charge transfer tendencies. 
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Elastic and electronic properties of fluorite R11O2 from first principle 



npy>i<Hi Ta e/ieKTpoHHi B/iacri/iBOCTi <|)/iK)opMTy Ru0 2 
3 nepuji/ix npi/iHi^nniB 

3.flx. flm 13 , A.M. rycP, TycP, flx. /P, 3. BainP, K. /liy^l P.ct>. /liHrxeP, X./J. flm 13 

npnpoflHH4nPi cfiaKy/ibTeT, TexHO/iormHufi yHiBepcuTeT MxeipflHy, 310023 XaH4>K0y, KHP 

tt>i3M4HMW cfiaKy^bTeT, Ka/ii(j)opHificbKnw flepxaBHi/iw yHiBepcuTeT, HopTpifl>K, Ka/iitj)opHm 91330-8268, CLUA 

<t>i3M4HMpi cfiaKy^bTeT, neflaroriHHwti yHiBepcuTeT HePiL^3flHK), HefiL(3flH 641112, KHP 

Ko/iefl>K MaTepia^03HaBCTBa i xiMiHHOi' mxeHepii, XaPiHaHCbKa L(eHTpa.nbHa periOHa/ibHa floaiiflHa 

yia6opaTopia 3 ni/iTaHb yTM/ii3ai4iT Si-Zr-Ti, XaPiHaHCbKnti yHiBepcuTeT, Xatkoy 570228, KHP 

tt>i3i/i4Hi/iPi cfiaKy^bTeT, MyHL4WHCbKwti TexHO^oriHHwR yHiBepcuTeT, Myt-mnH 400050, KHP 

tt>i3M4HMpi cfiaKy^bTeT, neflaroriHHWM yHiBepcuTeT ry^4>K0y, ryfiflH 550001, KHP 

iHCTHTyT aTOMHoT i MO/ieKy/wpHoT cf>i3MKH, Cn4yaHCbKnPi yHiBepcuTeT, HeHfly 610065, KHP 

npy>KHi, TepM0flWHaMi4Hi Ta eyieKTpi/i4Hi B/iacTHBOCTi cfi/iioopnTy RUO2 npi/i BucoKOMy TUCKy floaiiflxyroTbCfl 3a 
flonoMoroK) TeopiT tj)yHKL(ioHa.ny rycTHHH 3 n/iocKOXBH/ibOBWM nceBflonoTeHL(ia;iOM. OnTHMi30BaHi napaiweTpw 
rpaTKM, npyxHi cra/ii, o6'cmhwm MOfly/ib i MOfly/ib 3cyBy y3roflxyK)TbCfl 3 iHiunMW TeopeTH4HWMH 3HaneHH?iMH. 
ct>a30Bi/m nepexifl 3 MOflwefiiKOBaHoro cfi/iKiopnTy flo cfi/iroopnTy e npn 88 GPa (Ha6/iw>KeHHfl /lOKa/ibHOi rycTWHW, 
LDA), hh npn 1 15.5 GPa (y3ara/ibHeHe rpafliCHTHe Ha6/iw>KeHHfl, GGA). TaKOx floaiiflxeHO Mo^y/ib lOHra i Koecfii- 
L4i£HTM /laMe npw bwcokwx TUCKax. CTpyKTypa e CTa6i/ibHOK) p,nn tuckIb p,o 120 GPa, hklho 0641/icmoBa™ npyxm 
cra/ii. KpiM Toro, flooiiflxeHO TepMOflUHaMiHHi B/iacTHBOCTi, BK/iK)HaK)HH TeMnepaTypy fle6afl, TenyiocMHicrb, 
Koe<))ii4ieHT TenyiOBoro po3iunpeHHfl, napaMeTp TpK)Hati3eHa i KoetjjiLjieHT nyaccoHa. B e/ieKTpoHHiw CTpyKTypi 
c|)/iK)opnTy R11O2 3HaPifleHO Ma/iy 30HHy mi/inHy i LunpnHa 30hw 3pocTae i3 thckom. TaKox, npeflcraB/iem Me- 
xaHi^Hi Ta eyieKTpoHHi B/iacTHBOCTi fleMOHCTpyiOTb, mo npupofla 3B'H3yBaHHfl e KOM6iHai^icK) KOBa/ieHTHoro, 
iOHHOro i MeTa/iHHOro BioiafliB. 

KjiiOMOBi c/iOBa: nepwi npuHLf^nn, e/ieKTpoHHa crpyKTypa, npyxt-iicTb, TepMOfli/iHaMmHiCTb 
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